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ABSTRACT 

Combinatorial interactions of transcription modula- 
tors are critical to regulate cell-specific expression 
and to drive direct cell reprogramming (e.g. 
trans-differentiation). However, the identification of 
key transcription modulators from myriad of candi- 
date genes is laborious and time consuming. To 
rapidly identify key regulatory factors involved in 
direct cell reprogramming, we established a multi- 
plex single-cell screening system using a fibroblast- 
to-monocyte transition model. The system 
implements a single-cell 'shotgun-transduction' 
strategy followed by nested-single-cell-polymerase 
chain reaction (Nesc-PCR) gene expression analysis. 
To demonstrate this, we simultaneously transduced 
18 monocyte-enriched transcription modulators in 
fibroblasts followed by selection of single cells ex- 
pressing monocyte-specific CD14 and HLA-DR 
cell-surface markers from a heterogeneous popula- 
tion. Highly multiplex Nesc-PCR expression analysis 
revealed a variety of gene combinations with a 
significant enrichment of SPI1 (86/86) and a novel 
transcriptional modulator, HCLS1 (76/86), in the 
CD14 + /HLA-DR + single cells. We could further dem- 
onstrate the synergistic role of HCLS1 in regulating 
monocyte-specific gene expressions and phago- 
cytosis in dermal fibroblasts in the presence of 
SPI1. This study establishes a platform for a multi- 
plex single-cell screening of combinatorial tran- 
scription modulators to drive any direct cell 
reprogramming. 



INTRODUCTION 

A defined combination of transcription modulators 
is key to target-directed cell reprogramming or trans- 
differentiation. Retro/lentiviruses have enabled genetic 
approaches in cultured primary cells to ectopically 
express genes of interest (1-6). Viral screens using 
focused overexpression libraries can identify putative 
reprogramming factors, for example, by identifying a 
defined set of transcription factors that are selectively 
required to induce pluripotent stem cells (2-4,7) or 
trans-differentiate from mature somatic cell into another 
mature somatic cell (1,8—11). Despite the power of viral 
technologies, pooled screens for cell reprogramming have 
generally adopted a one-by-one elimination approach or 
testing for exhaustive combinations and interrogating 
phenotypes in cell culture from transgenic mice (1,3,10) 
and/or aberrantly cloned cell lines (12). These approaches 
require both considerable specialty and a substantial in- 
vestment in time. 

We therefore implemented a rapid and direct approach 
by allowing multiple viruses to infect at random and 
identify specific viral transcripts that give rise to 
target-specific expressions at a single cell level. Our focus 
was identifying essential interplay of transcription 
'modulating' factors (TMs) through lentivirus that were 
selectively expressed at a single-cell level. This type of 
screen holds promise for discovery of novel combination 
of TMs for regenerative medicine (1,13) or for any 
screening assay, which requires selectable phenotype 
such as biomarkers and cancer stem cells (14,15). 

We and others previously described that the genetic 
interactions is key to defining the transcription regulatory 
network of human monocytes (16,17). We thus hypothe- 
sized that introducing the transcriptional network of 
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monocytes into fibroblasts would promote trans- 
differentiation into functional monocytes characterized 
by the expression of both CD 14 and HLA-DR surface 
antigens. In this study, we transduced a pool of 18 
monocyte-enriched transcription modulators into human 
dermal fibroblasts and isolated CD 14 and HLA-DR- 
positive single cells after 2 weeks. A highly multiplex 
nested-single-cell-polymerase chain reaction (Nesc-PCR) 
analysis showed heterogeneous infection patterns but 
revealed significant enrichments of SPI1 and HCLS1 tran- 
scripts when selected for monocyte markers. On the basis 
of our single-cell screening method, we further demon- 
strate that SPI1 promoted and HCLS1 enhanced marker 
expressions and function specific to monocytes in dermal 
fibroblasts. 

MATERIALS AND METHODS 

Complementary DNA and virus preparation 

Gateway-compatible human full-length complementary 
DNA (cDNA) entry clones derived from RIKEN BRC 
clone bank [http://www.brc.riken.jp/ (23 July 2012, date 
last accessed)], Invitrogen (Carlsbad, CA, USA) and 
OpenBiosystems (Huntsville, AL, USA) were recombined 
into 150ng of pENTR lentivirus vector (CSII-EF- 
RfA-IRES2-VENUS; Supplementary Figure SI) over- 
night at room temperature using Gateway LR clonasell 
enzyme mix (Invitrogen). After 2h of Protease K treat- 
ment at 55°C, recombined plasmids were transformed 
into OneShot® Stbl3 competent Escherichia coli 
(Invitrogen) following manufacturer's protocol. Plasmids 
derived from single colony were expanded and purified 
using PureYield Plasmid Midiprep System (Promega). 
For every 8.5 ug of the plasmid, 5 ug of HIV-gp and 5 ug 
of VSV envelop genes were co-transfected onto 4x 10 
293T cells (prepared the day before at 37°C, 10% CO,) 
using FuGeneHD (Roche) in OPTI-MEM (WAKO) 
medium containing 5% FBS at 37°C 5% C0 2 . Between 
24 h and 72 h of incubation, supernatant-containing virus 
were collected and centrifuged at 19 400rpm for 2h at 
20°C. The pallet was then dissolved in 100 ul HBSS 
buffer (WAKO), followed by titer check and stored at 
— 80°C freezer for later use. 

Virus titer check 

One microliter of concentrated virus from the previous step 
was serially diluted (1:1000, 1:3000, 1:9000, 1:27 000 and 
1:36000) in MEM-a containing 10% FBS/L-glutamine/ 
antibiotics (WAKO) and transduced onto 3000 of 293 T 
cells seeded in a black clear-bottom 96-well plate (BD) 
for 3 days. The nucleus was stained using Hoechst 
(Invitrogen) for 30min at 37°C 5% C0 2 , and the plate 
was subjected to Cellomics ArrayScan XTi Reader 
(Thermo Scientific) for image analysis. FITC and Hoechst 
filters were used to detect Yellow Fluorescent Protein de- 
rivative (Venus) and nuclear staining, respectively, at 10 x 
magnification, and 10 images were taken per each well. 
Cellomics bio-application protocol 'SpotDetector' was 
used to segment Hoechst+ nucleus and to count Venus 



overlapping cells. The titer was calculated by 

Titer (c.f.u) = ([dilution factor x % of Venus+cells 
x initial cell number x volume per well]/ 
volume of concentrated virus). 

Cell culture and virus transduction 

Human CD14 + primary monocytes were purchased from 
Lonza, and THP1 cells (16) were cultured in monocyte 
medium containing RPMI1400 (WAKO) containing 
10% FBS/5M HEPES/2.5M sodium pyruvate/3.4 ug 
P-mercaptoethanol/L-glutamine/antibiotics (GIBCO). 
Human dermal fibroblasts (Passages 5-8) were purchased 
from RIKEN BRC and prepared at 5 x 10 4 cells in culture 
media: MEM-a containing 10% FBS/L-glutamine/antibi- 
otic (WAKO) at 37°C 5% C0 2 1 day before virus trans- 
duction. A total of 10 multiplicity of infection (MOI, 0.55 
MOI per TM) were pooled into 3 ml of culture media 
containing 80mg/ml of Polybrene (SIGMA) to increase 
transduction efficiency. The virus was transduced and 
incubated for 2 weeks in 37°C 5% C0 2 with MEM-a con- 
taining 10% FBS/L-glutamine/antibiotic in the period of 
first 5 days and in the monocyte medium for the remaining 
9 days. 

Immunostaining and cell sorting 

Two weeks post-virus transduction, cells were detached 
using Accutase and stained with human anti-CD 14/ 
AlexaFluora700 (BioLegend) and human anti-HLA-DR/ 
APC (BD) antibodies for 30min in ice. Fluorescence- 
activated cell sorting (FACS) analysis was performed on 
BD FACS Aria II (BD) for sorting Venus + /CD14 + / 
HLA-DR + single cells. To ensure a successful capture of 
single cells, we verified the cell-sorting accuracy before 
every experiment by observing the cells under a fluorescent 
microscope to confirm the alignment (center of the well) 
and to detect a single cell body. We detected 72 ± 10% 
positive single cells. Thereafter, single cells were sorted 
directly into 96-well plate containing the pre-amplification 
mixture. Pre-amplified samples positive for ACTB expres- 
sion were used for the subsequent analysis. 

Candidate gene extraction 

Illumina Human WG-6 version 3.0 microarrays were used 
to analyze differentially expressed transcripts between 
human dermal fibroblasts and CD 14+ monocytes (18) 
(GSE27304). After quantile normalization, top 20 
most differentially expressed transcripts annotated as 
'regulation of transcription 1 [based on Gene Ontology 
(19)] in monocytes when compared with fibroblasts were 
selected for lentivirus production (Student's r-test: 
P value < 0.001). The viral transcripts, GAS7 and 
EGR2, could not be detected using Nesc-PCR; hence, 
they were removed from the virus pool. 

Nested-single-cell PCR 

The primer mix for pre-amplification consisted of 20 uM 
of sense primers targeting ORF region of each gene (total 
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36) plus 100 uM of Attb2 anti-sense primer or 20 uM of 
anti-sense primer for endogenous transcripts for primary 
fibroblasts and monocytes. The primer mixture is further 
diluted to 200 nM and mixed with OneStep qRT-PCR 
2 x Reaction Mix, SuperScriptlll RT/Platinum TaqMix, 
RNAse OUT (Invitrogen) and DNAse/RNase-free dis- 
tilled water (GIBCO) to total volume of 9 ul per 
reaction. After single-cell sorting, the samples were 
incubated at 55°C for 25min, 95°C for 2min and 18 
cycles of 95°C for 15s and 60°C for 4min. The final 
samples were diluted 1:5 with DNase/RNase-free water 
and stored in — 20°C for later use. Primer mixture for 
PCR amplification consisted of 20 uM of both sense and 
antisense primer pairs targeting ORF region of each gene 
plus 0.3 uM of UPL probe for specificity (ROCHE). 
Additional 2x assay loading reagent (Fluidigm) and 
water were added to adjust final volume of 5 ul per 
reaction. Diluted pre-amplified cDNA (1.25 ul) is mixed 
with 2x FastStart Universal Probe Master (ROX; 
ROCHE) and 20 x GE Sample loading reagent 
(Fluidigm) and water to final volume of 5 ul per 
reaction. The primer mix and the sample mix were 
loaded onto 48.48 Dynamic Array (Fluidigm) and 
primed and mixed using the MX IFC controller 
(Fludigm). The dynamic array was further subjected to 
BioMark Array System (Fluidigm) at 95°C for lOmin 
and 40 cycles of 95°C for 15min and 60°C for 1 min. 

Single-cell data analysis and statistics 

Expression data from BioMark (Fluidigm) were processed 
in R Bioconductor (20). To illustrate gene expression and 
for statistical analyzes, we first set all Ct values >35 to 35 
and inversed the expression by subtracting the values by 
35. To test the statistical significance of marker enrich- 
ment analysis, we performed chi-square test to compare 
observed data (detected number of exogenous TM) with 
data we would expect to obtain (marker+ cells and 
marker — cells). The data are presented in terms of 
percentages. 

qRT-PCR 

Reverse transcription of total ribonucleic acid (RNA) was 
achieved with PrimeScript™ Reverse Transcriptase 
(Takara) and random hexamers in accordance with the 
manufacturer's protocol. Beta actin (ACTB) messenger 
RNA was used as a control for data normalization. The 
PCR primers used for each gene in this analysis are 
provided in Supplementary Table S2. PCR amplification 
was performed on an AB1 PRISM® 7500 Sequence 
Detection System (Applied Biosystems). For amplifica- 
tion, SYBR Premix Ex Taq™ II (Takara) was used as 
instructed in the manual. Changes of gene expression were 
determined using the 2— AACt method. 

Phagocytosis assay 

Combinations of vector control, SPI1 and HCLS1, were 
transduced onto human dermal fibroblasts seeded at 
5000 cells in a black clear bottom 96-well plate for a 
period of 20 days in RPMI1400 (WAKO) containing 
10% FBS/5 M HEPES/2.5 M sodium pyruvate/3.4 ug 



p-mercaptoethanol/L-glutamine/antibiotics (GIBCO) at 
37°C 5% C0 2 . Each plate contained six replicate wells 
of the same condition. The experiment was performed 
two independent times. One vial (2mg) of the pHrodo 
BioParticles fluorescent particles (Invitrogen) were 
diluted in 4 ml of uptake buffer consisting of HBSS 
(WAKO), 20 mM HEPES, pH 7.4 and replaced the 
culture medium with 100 ul of the prepared pHrodo 
BioParticles suspension. The plate was incubated in 
37°C 5% C0 2 for 7h, and cells were stained with 
Hoechst (Invitrogen) for lh at 37°C 5% C0 2 . The 
images were taken using the Leica XT microscope 
(Leica) and quantified using Cellomics ArrayScan XTi 
Reader (Thermo Scientific). 

siRNA transfection and RNA extraction 

Reverse transfection of 1 x 10 6 THP1 cells in each 60-mm 
cell culture dish was performed with 20 nM of each stealth 
negative control siRNA, SPI1 siRNA (Stealth RNAi: UG 
GUGCCCUAUGACACGGAUCUAUA) or HCLS1 
siRNA (Silencer® Select, Invitrogen; siRNA ID s6484) 
in Opti-MEM and Lipofectamine 2000 (Invitrogen), ac- 
cording to the manufacturer's instructions. Total RNA 
was extracted 48 h after transfection using miRNAeasy 
kit (Qiagen) according to manufacturer's instructions. 



RESULTS 

Single-cell screening system 

To rapidly screen for TMs involved in the direct monocyte 
specification, we designed an experimental workflow 
that integrates lentivirus technology (5,6), FACS and 
single-cell profiling analysis (21-23) (Figure 1). Human 
primary dermal fibroblasts were simultaneously trans- 
duced with a lentivirus pool of 18 monocyte-enriched 
TMs, selected based on gene-specific expression in 
primary CD 14+ monocytes when compared with dermal 
fibroblasts (Supplementary Table SI, see Section 
'Materials and Methods'). The stochastic nature of the 
virus infection allows for a mixed combination of TMs 
to be expressed per infected single cell, giving rise to 
target-specific expressions from a sub-population that 
accompanies the defined combination of monocyte- 
specifying TMs. A multiplex expression profiling of both 
exogenous and endogenous transcripts from single cells 
allows for the immediate identification of TMs essential 
for direct cell reprogramming. 

First, to ascertain accurate detection of exogenous 
genes, single cells expressing Venus fluorescent protein 
(i.e. 18-virus pool infected fibroblasts) were individually 
sorted using FACS and subjected to a multiplex Nesc- 
PCR for the detection of 18 virus-specific transcripts 
(Figure IB). Based on our statistical analysis, the 
single-cell transduction followed the Poisson distribution 
and the number of detected genes correlated with the half 
of estimated number of virus particles per cell (Figure 1C). 
This suggests that a minimum of two copies of viral tran- 
scripts can be detected, and 10 MOI of the pooled virus, or 
0.55 MOI per gene, was optimal to detect multiple target 
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Figure 1. Workflow of single-cell screening and Nested-single-cell-PCR (Nesc-PCR). (A) Lentviruses encoding multiple transcription 'modulating' 
factors (TMs) are pooled and transduced into human dermal fibroblasts for a period of 2 weeks. Because of the stochastic nature of virus infection, a 
variety of TM combinations are expressed at random per single cell. Cells expressing a defined set of TMs give rise to the expression of target-cell- 
specific markers (e.g. CD14 and HLA-DR) allowing them to be individually sorted using FACS. Subsequently, single cells are profiled using a 
Nesc-PCR gene expression analysis by reverse transcribing (RT) the exogenous transcripts with a virus-specific primer (primer 1) followed by 18 
cycles of pre-amplification with gene-specific forward primers (primer 2) and primer 1. The cDNA products are then subjected to a microfluidic 
qPCR for the identification of exogenous transcripts using gene-specific sense (primer 2) and anti-sense (primer 3) primer pairs and fluorescence 
probes (yellow star) (B). (C) The single cells transduced with pooled lentivirus (18 factors) at 5 MOI and 10 MOI were sorted and profiled using 
Nesc-PCR. At 10 MOI, >30% of single cells expressed 4 out of 18 TMs, while 20% of single cells expressed 2 out 18 TMs at 5 MOI. 



genes from single cells since >50% of the population con- 
tained average number of two to six genes. 

Single-cell gene expression analysis of virus-infected 
fibroblasts 

CD 14 and HLA-DR proteins have been well characterized 
to play a role in LPS response and antigen presentation in 
monocytes, respectively (24,25). To determine which of 
the 18 factors were critical for activating monocyte gene 
expression, we labeled the transduced cells using CD 14 
and HLA-DR antibodies. Two weeks post-infection at 
10 MOI, 18 ± 2% of the fibroblasts showed Venus fluor- 
escence (Figure 2A). Predictably from a pooled virus, a 
minute amount of Venus + cells were positive for both 
CD14 and HLA-DR (0.17 ± 0.1%; Figure 2B) but was 
sufficient to obtain single cells for transcript profiling 
using Nesc-PCR (M+, n = 86). To analyze the enrichment 
of TMs in marker-positive single cells, we additionally 
sorted an equal number of single cells from Venus + popu- 
lation as a baseline control (Figure 2A, M— ). 
Additionally, fibroblasts transduced with vector control 
(10 MOI) failed to express CD14 and HLA-DR surface 
markers (data not shown). 

Individually sorted single cells were profiled specifically 
targeting 18 viral-TM genes, 16 endogenous monocyte 



markers and 1 internal control (Supplementary 
Table S2). The final 35 monocyte-gene panel was used 
to measure by means of Nesc-PCR. To further address 
the effect of exogenous factors on fibroblasts, we addition- 
ally isolated and analyzed single cells from 'unperturbed' 
fibroblasts and human CD14 + primary monocytes (« = 40 
each, Supplementary Figure S2A). These individual cells 
were profiled using the same 35-gene panel but targeting 
18 endogenous TMs. 

We then analyzed the ability of Nesc-PCR gene expres- 
sion analysis to distinguish different cell populations 
based on exogenous TMs and endogenous marker expres- 
sions. Figure 3A reveals a cluster of 86 CD14 + /HLA-DR + 
cells (M+) and 86 Venus + cells (M-). Most notably, single 
cells positive for CD 14 and HLA-DR markers detected 
SPI1 (86/86, 100%) and 76 out of 86 expressed HCLS1 
(88%), which were statistically significant when compared 
with the M— population where only 8/86 and 10/86 single 
cells were positive for SPI1 and HCLS1, respectively 
(chi-square test, P value < 0.001 for both). Interestingly, 
though eight out of nine single cells expressing exogenous 
SPI1 in the M— population expressed CD14 transcripts, 
and two out of nine expressed HLA-DRA transcripts, but 
no detection of HLA-DRB transcript was observed. 
Although the sample size was small, the analysis 
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Figure 2. Flow cytometry sorting of CD14 and HLA-DR-positive single cells. Two weeks post transduction of 18 monocyte-TMs, 18% of fibroblasts 
expressed Venus fluorescent protein (A, 'M— '), out of which 0.17% was positive for both CD14 and HLA-DR markers (B, 'M+'). Eighty-six single 
cells from M+ and M— boxes were sorted for further analysis. 



suggests that SPI1 alone could induce the expression 
of CD14. In fact, when SPI1 alone was transduced into 
fibroblasts in an independent experiment, 12 ± 2% of 
SPI1 + cells expressed the CD 14 surface marker suggesting 
that SPI1 could elicit the expression of this marker but at a 
low efficiency (Supplementary Figure S3). 

In a drastic contrast to the M— population, the analysis 
of M+ population noticeably revealed that this subset, 
although transcriptionally heterogeneous, was almost ex- 
clusively composed of cells expressing high levels of genes 
characteristic of mature monocytes (Figure 3A). Single 
cells selected for CD14 and HLA-DR surface proteins 
revealed 100% detection of CD14, HLA-DRA and 
HLA-DRB marker transcripts, as expected. But add- 
itional markers such as CSF1R (85%), ILT4 (73%), 
TYROBP (63%) and CD 163 (52%) transcripts were 
highly enriched in the M+ population when compared 
with the M— population (chi-square test, P value < 0.05) 
and significantly induced when compared with the 
basal expressions in dermal fibroblasts (Supplementary 
Figure S4). 

SPI1 and HCLS1 as key determinants of monocyte 
expression 

To further investigate the key regulatory determinants of 
monocyte specification, we performed eigenvector-based 
analyses between M+, M— , fibroblasts and monocytes 
to visualize the heterogeneity in single cells, as reflected 
by the 2-axis of the principal component analysis (PCA; 
Figure 3B). The analysis found that the first principal 
component (explaining 27.57% of the variation) grouped 
the primary monocytes and the M+ cells but separated 
fibroblasts and M— cells. These distinct clusters indicate 
that M+ cells deviated from the original fibroblast state 
and transitioned into a state comparable to that of CD14 + 
monocytes. However, a closer examination of the second 
component narrowly distinguished the monocytes from 
the M+ cells, indicative of partial (or incomplete) induc- 
tion of monocyte network in fibroblasts. The PCA also 
showed a wide scatter of M— single cells largely due to 
random virus transductions. However, a small number of 



single cells clustered closely to the M+ population largely 
contributed by positive detection of SPI1 and CD 14 ex- 
pression but failed to be selected by FACS due to lack of 
HLA-DR protein expression. 

To identify the causal genes that contribute to the sep- 
aration of each cell population, we performed the corres- 
pondence analysis (Co A; Figure 4C). Complementary to 
PCA, the CoA identified SPI1 and HCLS1 as major de- 
terminants to discriminate the monocytes from the 
original fibroblasts. Although the primary monocytes ex- 
pressed SPI1 and HCLS1, the exogenous expression in the 
M+ population was 6-fold higher (in the case of SPI1), 
contributing to the separation from the primary mono- 
cytes on the second axis. Interestingly, the single-cell 
expression profile of primary monocytes revealed that 
27% and 22% monocytes expressed CEBPA and IRF8, 
respectively, which is in contrast to the ubiquitous expres- 
sion of SPI1 and HCLS1 in all 40 monocyte single cells 
profiled (Supplementary Figure S2B). The heterogeneous 
expressions of CEBPA and IRF8 transcripts may suggest 
a dynamic balance of expression to maintain monocyte 
phenotype but are not necessary required to promote 
key monocyte-specific markers. 

SPI1 is required and HCLS1 enhances expression and 
function specific to monocytes 

To further explore the effects of combinatorial TMs, we 
ectopically expressed vector control, SPI1, SPI1+HCLS1 
in fibroblast and measured the differential expression 
using the conventional qRT-PCR method. CSF1R is a 
well-studied monocyte membrane protein acting as a 
receptor for colony stimulating factor 1 , and deregulation 
of this gene is a hallmark of many tumors (26). Although 
SPI1 alone could induce the expression of this monocyte 
marker, the signal intensity was greatly enhanced when 
fibroblasts were co-infected with HCLS1 (Figure 4A). 
Similar expression patterns were also observed for 
TYROBP and CD163. However, HCLS1 could not 
enhance the transcript level of CD14. When the induced 
expressions were compared with the basal expression of 
monocytes, the cell transformation was only partial 
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cells transduced only with vector control or HCLS1. Scale bar = 100 urn. (D) Gene-specific knockdown (KD) of SPI1 and HCLS1 in THP1 cells 
revealed significant suppression of IRF8 and CEBPA transcripts, respectively. Knockdown of SPI1 also suppressed CSF1R, but induced CD 163 and 
CD14 expression, which was also observed when HCLS1 transcript was ablated. P value (Student's /-test) *< 0.001, when compared with vector 
control. 



(Figure 4B). This incomplete observation was also seen in 
the single-cell PCA analysis (Figure 3B), indicating that 
the conversion is heterogeneous and partial, although 
some single cells clustered closer to the primary monocytes 
than others. Concurrently, it is also plausible that only a 
subpopulation of the cells undergoes cell transformation, 
diluting the expression levels and thus masking the 
induced expression. 

Next to test whether SPI1 in conjunction with HCLS1 
could also enhance the functional role of monocytes, we 
added pHrodo™ Red E. coli BioParticles® conjugate 
beads, which is pH-sensitive Rhodamine dye to fibroblasts 
transduced with lentivirus vector control, SPI1, HCLS1 or 



both (Figure 4C; Supplementary Figure S5). Since the 
fluorescent dye of the beads specifically expresses in 
phagosomes, we could reliably detect the ingested beads 
using the fluorescent microscopy without trypsinization of 
the cells. We observed a clear detection of Rhodamine- 
positive cells after transducing SPI1 and the SPI1/ 
HCLS1 combination. The image analysis further 
revealed that cells transduced in combination of SPI1 
and HCLS1 ingested 33% more beads when compared 
with SPI1 alone (Supplementary Figure S5B). 

To better understand the causal effects of SPI1 and 
HCSL1 involved in monocyte specification, we targeted 
the transcripts using siRNA oligos in THP1 (human 
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acute monocytic leukemia cell line). THP1 has been 
reported to be a suitable model to study monocyte 
network (16) and also has been demonstrated to yield 
high efficiencies in siRNA transfections (16,27). Since 
IRF8 and CEBPA have been implicated as key regulators 
of monocyte network, we specifically measured the expres- 
sion changes after the knockdown (16,18,28,29). 
Interestingly, ablating SPI1 and HCLS1 transcripts 
in THP1 cells could significantly repress IRF8 and 
CEBPA expression by 2.71-fold and 2.56-fold, respectively 
(Figure 4D), as well as CSF1R by 6.41 -fold — in the case of 
SPI1 (Figure 4E). Surprisingly, however, THP1 deficient 
in SPI1 augmented CD 163 and CD 14 expression, and a 
similar expression pattern was also observed for the 
HCLS1 knockdown. Although the exact mechanisms is 
unclear, ablation of SPI1 and/or HCLS1 in the established 
monocyte network may not be sufficient to suppress 
marker genes such as TYROBP ('gene redundancy') or 
may induce gene expressions to maintain network homeo- 
stasis ('network robustness'). Taken together, our results 
revealed that SPI1 is minimally required to promote, in 
part, monocytic expression and phagocytic function in 
human dermal fibroblasts but can be greatly enhanced 
through the co-expression of HCLS1. 



DISCUSSION 

Here we described a rapid and direct method that facili- 
tates the use of lentivirus and Nesc-PCR for identifying 
key genes involved in trans-differentiation. By relying on 
the random nature of virus infection, we implemented a 
single-cell 'shotgun-transduction' system that enabled the 
prompt detection, isolation and identification of combina- 
torial genes that are required to transform fibroblasts into 
a desired target cell type. 

By taking advantage of viral transcripts and the 
single-cell detection method, we noted SPI1 and HCLS1 
as a defined set of transcription modulators responsible to 
switch from human dermal fibroblasts to functional 
monocyte-like cells. SPI1, an important regulator of 
myeloid and B-lymphoid cell development, also acts as a 
key modulator of cell transformation toward monocytes 
together with CEBPA (9). The detection of SPI1 by the 
single-cell screening system decisively validates our 
method in identifying the strongest modulator in the 
group. However, we could additionally demonstrate that 
CEBPA was dispensable in that SPI1 alone could induce 
the expression of CD 14 and HLA-DRA and promote 
phagocytosis in human fibroblasts. Although HCLS1 
alone could not induce monocyte trans-differentiation, 
the protein acted as a strong modulator of monocyte/ 
macrophage phenotype. HCLS1 has been implicated in 
response to Cortisol and LPS, activating numerous func- 
tions of macrophages (30,31). Although the exact mech- 
anism is unknown, our data and the published reports 
strongly suggest that HCLS1 expression correlates posi- 
tively with enhanced function of monocytes/macrophages 
possibly via regulating CEBPA expression. 

In this study, we have validated a highly scalable 
approach for screening lentiviral libraries for direct cell 



reprogramming. We used cell-specific markers to isolate 
single cells; however, this method can virtually be applied 
to any characteristics that allow separation of phenotyp- 
ically distinct cells, such as cell proliferation dyes or func- 
tional biomarkers for the identification of transcription 
modulators. Overall, this study demonstrates a rapid 
single-cell screening workflow to readily identify and char- 
acterize new combinatorial factors, applicable to induce 
various cell-reprogramming pathways. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online: 
Supplementary Tables 1 and 2 and Supplementary 
Figures 1-5. 
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